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Abstract —Nowadays technological developments, the use 
of lasers in production is increasing and plays an 
important role due to low cost and high accuracy. The heat 
transfer, over the course of laser-metal interplay, has a 
great importance in metal forming. In this study, different 
types of materials were investigated in order to designate 
the temperature distributions inside material and on the 
material surface versus the thermodynamic properties of 
the material used and then the temperature distributions 
obtained from the analysis were compared each other. 

In addition, the heat transfer is occurring during the 
interaction of the laser power of 1.10 10 W/m 2 and 5.10 10 
W/m 2 with laser power intensity in two main groups using 
different materials these are steel, nickel, tantalum and 
titanium, and numerical results are obtained using the 
finite-difference method. In the first step, a solution is 
obtained by electron kinetic theory according to the basic 
heat transfer. In the second step, since heat convection is 
formed after material has reached the melting point. Using 
electron kinetic theory model for convection solutions have 
been obtained. Moreover, the temperature distribution that 
occurs during the laser metal interaction was studied by 
variation of the time chart and the material depth. As a 
result of the study, material's surface at the correct 
temperature of liquid phase change material andincreased 
depth in the direction perpendicular to the electro-kinetic 
theory approach is further demonstrated by the decrease 
in the first manner and then remains constant in 
exponential phase change temperature. In addition to this 
the analysis results, the substrate temperature increases, 
the change in phase in the material becomes smaller and 
smaller. 

Keywords — Kinetic Theory, Pulsed Laser-Metal , laser. 

I. INTRODUCTION 

A laser is a device electromagnetic stimulated 
emission of radiation that continues to happen for a while 


instantly optical amplification of light emission. The term 
"laser" with stimulated emission of radiation "has emerged 
as an acronym for light amplification [1]. In today's 
developing technology,because of the delicate processing 
quality of the laser material, laser has many uses such as 
engineering, industrial processing [2], medicine, scientific 
research [3], electronic and communication field [4]. The 
basic principle of laser devices is the output photon energy 
from the light source by passing it through a medium, of 
electrons from atoms very and available beams are very 
different motion in a single direction is to obtain a new 
beam [5]. The process describing the laser machining 
necessitate a knowledge of physical properties of the 
material and the machining ambient, such as heat capacity, 
thermal diffusivity, gases and pressure in the processing 

environment [6]. In the laser machining -process, 

various phenomena that occur when the electromagnetic 
radiation that comes out on the surface of a material come 
into the form of scattering, absorption and transmission. 
The absorption of radiation leads to various effects in 
materials, such as evaporation, possible plasma formation, 
molten layer formation, etc., which constitute the pressuies 
of laser material processing techniques [7]. 

Laser metal interacts with the free electrons in the 
transmission band to the surface of the photons that hit the 
surface of the material. As a result, electrons with 
increasing kinetic energy collide with excited cage 
molecules at higher energies. The higher vibrational 
energy of the molecules leads to higher temperatures. 
These interactions result in phonon-phonon confusion as to 
provide thermal equilibrium All of these events give the 
definition of the heat transfer mechanism [8]. 

In the literature, there are many various studies that help to 
distribute the temperature due to laser material interaction 
in solids and model it using Fourier equation. Koc at al. 
1995, investigated the phenomena of heat transfer 
occurring for the laser processing power of 1.10 10 W/m 2 
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and 5.10 10 W/m 2 . At the same time heat conduction and 
convection were evaluated electron kinetic theory, in 
addition to this different types of materials analysed in 
order to determine the temperature distraction [9] Yilbas at 
al. 2006, analysed and developed the electron kinetic 
theory formulation for laser ultra-short pulse heating by 
estimating the comparatively low cage temperatures. 
Consequently, the developed formula is the analogous with 
found for improved formulation hyperbolic heating model 
[6]. Mansooretal. 2015, The effect on the laser short-pulse 
temperature was investigated for an aluminium thin film 
and the energy transfer was calculated using the Boltzmann 
equation in the film When they calculate the time- 
dependent change in equilibrium temperature, the laser 
pulse intensity in the subsystem continues. However, the 
equilibrium temperature is different for the electron and 
lattice subsystems associated with the phonon [10]. Koc et 
al. 2004, examining the 3-D laser heating model on four 
different materials (tantalum, iron, titanium, nickel) using 
the electron kinetic theory approach. For these materials 
the temperature distribution depends on the time. 
According to this information the surface temperature of 
the material is calculated and analysed. Asa result of the 
study, He found the temperature rises up to the point of 
melting at the material origin and by the reason of to heat 
transfer by conduction. When heating for these materials, 
the temperature distribution is reduced in the x direction 
[11]. Shuja and Yilbas at al. 2013, investigated and 
simulated laser which is induced melting of coated carbon 
steel which is kept constant at 7.5 mm and is contained with 
alumina, boron, tungsten carbide. As a consequence of this, 
surface temperature raised during the first heating phase of 
the laser pulse, However the temperature of cooling cycle 
is precipitous drop [12]. They studied laser metal 
interactions and focuses on process of development for 
flexible applications. In addition to this, investigated new 
approaches with laser metal interaction which is flexible 
tool for transparent electrodes of metal or synthesis of 
nanomaterials [1]. 

II. MATHEMATICAL MODELLING OF 
HFAT TRANSFER 

Heat transfer analysis is required to examine the 
temperature distribution during the laser metal interaction. 
Fourier's equations of conductivity used in heat transfer 
processes are insufficient during high-power laser-metal 
interaction. For this reason, there is a need for kinetic 
theory to study the heat transfer phenomenon during laser- 
material interaction more accurately [11]. 

This model can be explained as follows. Electromagnetic 
radiation applied from the outside is absorbed by free 
electrons in the transmission line of metals. These 
electrons, rise in energy, collide with atoms and leave some 


of their energy in atoms. Because of this mutual effect, the 
phonon energies of the atoms increase and cause the 
phonon energies to soar by forcing them to vibrate around 
the atoms around them. Some acceptances have been made 
to facilitate the resolution of the problem [9] [13]. 

■ It is assumed that dissipation of energy due to 
thermionic emission do not occur. 

■ Local equilibrium,^ presumed at material t = 0 
uniform temperature. 

■ It is accepted that the plasma has a continuous 
regime, the absence of nucleation and liquid metal 
flushing, the neglect of the physical change of the 
laser beam due to the laser optical system and 
plasma. 

■ It is supposed that electrons will not be affected 
by temperature changes of the mean free path 
[14]. 

According to kinetic theory, it is necessary to show regard 
to the free electron model in the quantum field in order to 
study the conduction of the heat. Electrons that move freely 
within the surface at a certain speed are called free 
electrons. Within a volume confined to a particular surface, 
electrons can change their energy according to the law of 
dispersion of the fume. The electrons are placed in the 
molecular cage and have phonon energy according to their 
vibrational motion [15]. When considering the lattice 
vibration, the displacement of an atom at the "s" point in 
the lattice is given by the following expression [9]. 

Ag e i(gs-cogt) 

Where; motion's amplitude is A g . At the T temperature, 


the approximate coupling energy related to this vibration 
mode is given by the Planck formula [9] [14]. The energy 
of the lattice vibration contingent upon the phonon number 
(n) and frequency (co). If the number of phonons (n) in the 
T temperature and in the thermal equilibrium is written 
depending on the Planck distribution; 

i 

<11>= e *<B/k B T _ 1 


or; 


<n> + - = -ctnh. 
2 2 


ha* 

2k^T 


Where; < > : indicates the average of the thermal balance 
kB : Boltzman canstant. 

<n> : The average number of photons in thermal 
equilibrium 


Nn+l / N n — 


g/zco/ T 


T =kB.T 
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It is written in the form. In this case, the ratio of "nth" to 
the total number of oscillations in the quantum state; 


N r 


7X 

s=o 


g-n/zoo/ t 

00 

y e -s>/zcfl/T 

s=o 


If the regulation of the mean excitation quantumnumber of 
the oscillation is made and x = exp (-/zcox), the following 
takes the form; 


Is.e s 



X 

(1-x) 2 


In the classic sense, each oscillation is considered to have 
“ <n>/zco « kBT ” energy. The average energy of an 
oscillation at frequency co is <h> /zco . The average energy 
for one atom vibration; 

<E> = h. cog / [exp (h.cog/k.T-1)] 

C v is Heat capacity of the oscillations which is written in 
the form. 


c„ = 


f dE} 
\dT j 


= N.k B . 


Ti.co 


fico! T 


\ T ) ' [e- ncol% -l ) 2 


Equation 3N is written for three-dimensional N atom 
oscillations. Because each atom has three degrees of 
freedom N atomic energy in a single direction due to 
tremble vibrations; 

E x — N x .kg.T 

3N equation can be identified as the phonon energy. For 
the most part, in a homogeneous material, thermal 
conductivity factor is described by assuming that heat 
conduction is continuous and the thermal conductivity 
factor in the single direction is given by the following 
relation [16] [17]. 


IX N.V x .k B .A, 

K= 6 

In this analysis, mathematical explanations of electrons 
moving in the material and reflected back from the 
surface are strong. This situation can be removed from the 
surface by imaginary mirror placed on the surface. The 
expected arrival of an electron at x distance without 
collision is exp (x/X) for x«A. The probability of 
collision at the end of an electron dx interval is dx/X 
According to the figure 1, the possibility of collision last 
for II, III, I ways of the electrons at dx distance is as 


follows. 


N/6 




a) Electron flow b) Mirror Acceptance 

Fig.l: Direction of movement of electrons in relation to energy transport 


The probability of total collisions of electrons in the surface section along the medium might be written as follows [8] [13]. 

x — s\ ds dx 
X XX 

Where the distance of the “x” atoms and “A” is the electron free path [18]. The lower limit of the integral (-oo) in this equation 
shows the electrons returning from the surface into the material. 


00 
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—00 
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Electron movement is described in figure 1. Electron phonon collision result, the net energy transfer from electrons to phonons 
is written as follows. 


00 

AE x,t= 1 


exp - 


■CO 


ds dx 

T'T 


E s, t Ex, t 


The Es,t and Ex,t sequences in this form are the energies of electrons and atoms at one point, respectively. The parameter f is 
the percentage of energy that the electron gives to the atom during the atomic collision and f = 10 -4 is taken [19] [20]. The 
change in the passage of heat from a homogeneous medium which is a formula of distance x is written by the following 
equation. 

dx 

Absorption coefficient is denoted by 5. The negative sign here states a decrease in beam intensity thanks to absorption of the 
positive magnitude [21]. If we integrate this equation; 

I = I G .exp(-.x) 

Here; Io is the maximum intensity of the incoming laser beam Because of this reason, using the afore function for the high 
level intensity of the incoming laser beam, the energy absorbed during dt time domain dx might be given as: 

AE X t | = l 0 -8.e _ ^' x . A.dt.dx 
abs 

The overall energy augment in metal during dt time at dt [22]; 


N - A - ( E X, t + dt - E xA dx = AE x, t + AE x, t 


abs 


For this reason, the total energy transfer equation is [23]: 


N.A.(E X t + dt _E x t)- dx = lo- 5 - ex P( - 8.x). A.dt.dx 


00 — 
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•M E s,t- E x,t)- dsAdtdx 


The "V" in this form refers to the electron velocity. If the energy transfer in electron atomic collision is written in terms of 
electron and atomic temperatures, the formula is arranged as follows [24]: 
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III. MATHEMATIC MODEL OF HEAT CONVECTION WITH KINETIC THEORY APPROACH 

If an atom moving outward from the surface is in the vapour phase, the applied counter forces slow down the atom If the 
atom’s initial kinetic energy allows escape, the atom is vaporized [25]. The escape will be equal to the latent heat of evaporation. 
The initial minimum value required for escape Vx [26]: 

^ m -V x 2 =U 0 


The velocity distribution of the molecules according to Maxwell law is written as follows: 

m 


f(V x ).dV x = 


2 . 71 . k B .T 


" 1/2 

-m.V x 2 

.exp 

_2.kB.T_ 


■dVy 


Hence, the number of molecules separating the surface during time dt is: 

N. f (Vx).d Vx. Vx.dt 


Here N is the number of molecules on the surface of the unit volume. All these atoms with velocities greater than Vmin will 
be able to move away from the surface. G is the unit field, and the unit is the number of atoms evaporating is written as follows: 

oo 


G = N 


. I f(x).V x 

V . 
min 

,dV x 




m 

1/2 ^ 

I.XP 

Arr >\/ 2 

-m.V x 2 

d(mV x 2 ) 

2. 71 .kg.T 

' m' 

_2.k B .T_ 

2 


G = N. 


r k B .T 1 

1/2 

u 0 “ 

I 2 . 7 i.m_ 

exp 

L k B .Tj 


1 2 

Where using this formula Vq = ~. ITI.V , Vs is velocity of the evaporation of surface and can be written as: 


G 

N 


V c 


k B -T 

2 . 7 t.m. 


1/2 


exp 


U 


o 


k B -T 


There is no restriction to prevent molecules from escaping inside the cage. However, there will be a gap in the cage area behind 
the unreleased molecules. This space is filled by the free energy and other free molecules. The molecules released in the cage 
immediately form a new equilibrium position. In this way, the molecule is absorbed by a void of metal and can be moved at 
the lattice level. The average velocity of moving free molecules in the cage is much lower than the thermal velocities Vth [8] 
[11]. The Vf rate of the free molecules is approximately; 


V f = V th .exp 


U 


k B .T 


Here U is the energy essential to obtain the free molecule. This is important to take into account the energy delivered by free 
molecules and holes. The number of free molecules in any unit volume (N f ) in the T temperature; 

u, " 


N f = N 0 .exp 


k B -T 


The amount of the voids contained in a unit volume [11]; 
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N h = N 0 .exp 


U, 


k B -T. 


Here U2 is the activation energy. In this case, N 0 is the number of cage places in the unit cell. As a consequence, 


Nf=N h 


But this acceptance is not exactly true. 


N + N = N 


The number of conjugated molecule shown in N bm . If p is the average free path for a molecule that collides with alternative 
molecules, there is a possibility of a collision at dx is [6]: 

dx 

|Ll J ft 


, , dx 

1-exp -— 


The possibility of collision in a space P. ; 


P h = 


N 


h 


N 


h 


exp 




k R .T 


Nbm+Nfm+Nh No+Nf m 


The probability of collision with a free molecule is P m ; 

Nfm Nf m 


1 + exp 


exp - 


Ui 


k B -T 


Pm = 


Ui 


k B -T 


Nbm + Nfm + Nh Nq + Nf m 


1 + exp - 


Ul 

k B .T 


The probability of collision with a compound molecule is P bm ; 

u 2 


1-exp 


p bm = 


k B .T 


1 + exp - 


Ui 


k R .T 


There are two types of energy carriers that transport molecules. These are free molecules and holes. If the number of free 
molecules that collide with a space is equal to the number of moving molecules that are separated from the cage by any period 
at the time dt. This expression describes the results of the transmission energy from any dx cross section due to the escaping 
mobile molecules of the appropriate number for the free molecule hole collision operation [27] In the elastic collision process, 
the re-captured molecule will transfer some energy. Transfer energy [24] [28] ; 

U(s) + 3.k B .<|>(s,t) 

Here; § (s, t) is the lattice temperature in the s position at time t. The amount of energy that escaped molecules will take with 
it 

Ux) + 3.k B .<|>(x,t) 

The energy balance in dx results amount of energy. 

Us)-Ux) + 3.k B .[ <t>(s,t) - <t>(x,t) ] 

Where U (s) is the activation energy in the free position of the molecule s and is usually a function of the temperature. But 


www.iiaers.com 


Page | 73 



































International Journal of Advanced Engineering Research and Science (IJAERS) 

httos://dx.doi.ora/10.22161/iiaers.5.10.9 


[Vol-5, Issue-10, Oct-2018] 
ISSN: 2349-6495(P) / 2456-1908(0) 


since the surface is surrounded by molecules with a small amount, there is a decrease in the surface Vo. The number of free 
molecules separated from ds to dx in the time dt after a collision, where the molecules held in the eye are transferred from the 
energy ds section to the dx section; 


Nsx- A. Uq X . dt. 


sx- 


ds 


Here; Usx is the relative velocity of the free molecules from s to x The probability of collision of these molecules first in the 
dx section; 

ds dx 


Nqy • A. Ucy. dt. ■ 


p p 


The number of collisions with the free molecules is: 


N SX - A. Uqy. dt. 


ds dx 


'SX- 


P p 

The number of the collisions with escaping molecules is [11]: 

. ds dx 


.exp 


The number of collision with the holes is: 
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The amount of energy that the molecule transfers during the collision with a vacancy is; 

U(s)-U(x) + 3.k B .[ 4>(s, t) — <j)(x, t) ] 

The amount of the energy transferred by the molecule during the collision with a bound molecule is: 

^.g.k B .[ <j>(s,t) — 4>(x,t) ] 


Herewith, the energy which is transferred to the cross section dxfrom all the fields in ds is: 

^-.U sx .h.k B .P fm .exp 
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[ Us)-lXx) + 3.k B .[ 4>(s,t)-<j>(x,t) ] ]ds 


Thus, taking into account the transmis sion equation involving convection effects, the transmis sion and convection equation fo r 
the kinetic theory that gives the temperature distribution is written as follows [8] [29] 

K ' f T(x,t) 
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Mathematical Modelling of Pulse Heating 

The output of the beam from the solid lasers peaks in a few millimetres. The series of these peaks rise and eventually result in 
a fall. In the pulsed heating process, p and £, are chosen appropriately. These two exponential functional analytical equations 
are used and parameters give rise time. II is the maximum amplitude determined. 

l 0 =h-[exp(-pt)-exp(-Ct)] 


In order to see the effect of the impulse in the temperature distributions, the impulse N is modified to remain equal to the total 
power density and is shown in Figure 2. 
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12 


14 


16 


Fig.2: The change in power density of the pulsed laser over time. 

Pulse Heating Analysed and Numerical Solution 

During the laser-material interaction, the distribution of temperature in the material was analysed in one dimension depending 
on the time. In this study, the analysis of the temperature distributions during the interaction of the laser power of 1.10 10 W/m 2 
and 5.10 10 W/m 2 with the laser power intensity of Steel, Nickel, Tantalum and Titanium materials was made according to 
kinetic theory. 

The temperature variations at the material surface during the interaction of steel, nickel tantalum and titanium materials, 
respectively, with laser intensities of power densities of 1.10 10 W/m 2 and 5.10 10 W/m 2 , found according to the electro kinetic 
theory approach which is given in Figures 3 and 4. 


Times of melting temperature of materials; 

1.1Q 10 W/m 2 Power Intensity 


5.1Q 10 W/m 2 Power Intensity 


Steel 0.014 psn, 

Nickel 0.0174 jusn, 

Tantalum 0.0224 psn 


0.0028 jLisn, 
0.0035 jLisn, 
0.00448 jusn 


Titanium 


0.0149 jusn 


0.00298 jusn 



Steel Nickel 

Fig.3: Time-dependent variation of the temperature distribution of the steel and nickel material surface during the 
interaction with the laser beam at the power densities ofl.10 10 W/m 2 and 5.10 10 W/m 2 , respectively, by electrokinetic theory 

approach. 
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ZAMAN 1*E-6 (sn) ZAMAN 1*E-6 (sn) 


Tantalum Titanium 

Fig.4: Time-dependent variation of the temperature distribution of the tantalum and titanium material surface during the 
interaction with the laser beam at the power densities of 1.10 10 W/m 2 and 5 JO 10 W/m 2 , respectively, by electro kinetic theory 

approach. 

As can be seen from these figures, at a power density of 1.10 10 W/m 2 , the temperature of the material surface step-up rapidly 
until the melting temperature, then the temperature proceeds constant until the phase change is completed at the melting 
temperature. At a power density of 5.10 10 W/m?, the temperature increases slowly until it reaches the melting temperature, then 
increases slowly after reaching the melting temperature. The reason for this is that in high power density phase change occurs 
in a shorter period. 



Fig.5: Variation of the temperature distributions in the material according to material depth found by the electro kinetic 
theory approach during interaction with the laser beam at the power density ofl.10 10 W/m 2 the steel and nickel. 



Fig.6: Variation of the temperature distributions in the material according to material depth found by the electro kinetic 
theory approach during interaction with the laser beam at the power density ofl.10 10 W/m 2 the tantalum and titanium. 
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Figures 5 and 6 show the temperature distributions of the material depth at various times according to the electro kinetic theory 
approach of the interaction of steel, nickel, tantalum and titanium materials with the power density of 1.10 10 W/m 2 power 
density, 5.10 10 W/m 2 power density, respectively, in figures 7 and 8. 




Fig.7: Variation of the temperature distributions in the material according to material depth found by the electro kinetic 
theory approach during interaction with the laser beam at the power density of 5 JO 10 W/m 2 the steel and nickel. 



Fig.8: Variation of the temperature distributions in the material according to material depth found by the electro kinetic 
theory approach during interaction with the laser beam at the power density of 5.10 10 W/m 2 the tantalum and titanium. 


As shown, material is coming reduction occurs at the 
correct temperature to the depth from the surface 
temperature distribution in the liquid phase at least one 
exponential by showing decreased and remained stable at 
the phase change temperature, then showed a rapid decline 
in exponential solid phase then continues as 
asymptotically. Phase change material in the surface 
temperature increases, the temperature remains constant, 
decreases gradually in the liquid layer. 

IV. CONCLUSION 

When the surface of the material upon temperature increase 
take into account, the results of analysis made by electro 
kinetic theory approach can be considered in three regions. 
The solid phase found in the first zone. Although it appears 
linear graph very slight increase in temperature continues 
until the melting point of the parabolic. The energy density 
at low temperatures of the liquid phase of the second, first 
remains constant and then increases gradually. In this case 
the energy of the electrons and the kinetic energy increases 


pausing of atoms that is defined and explained by the 
increase in potential energy is converted into latent heat 
needed to phase changes. However, when high 
concentrations of the phase change material in contact with 
the laser beam is formed within a shorter time. 

When analysed as a time-dependent increase in surface 
temperature; heat transfer coefficient and absorption 
coefficient of the small increase in surface temperature of 
the material was determined to be fas ter than large ones. In 
contrast, materials with large heat transfer coefficient, the 
surface temperature increase becomes slow and material 
temperature increases towards depth is more rapid. 
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